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GENERALIZED SIMULATION TECHNIQUE FOR TURBOJET 
ENGINE SYSTEM ANALYSIS 

by Kurt Seldner, James R. Mihaloew, and Ronald J. Blaha 
Lewis Research Center 

SUMMARY 

A nonlinear analog simulation of a turbojet engine was developed. The purpose of 
the study was to establish simulation techniques applicable to propulsion system dynam- 
ics and controls research. 

A schematic model was derived from a physical description of a J85-13 turbojet 
engine. Basic conservation equations were applied to each component along with their 
individual performance characteristics to derive a mathematical representation. The 
simulation was mechanized on an analog computer. 

The simulation was verified in both steady- state and dynamic modes by comparing 
analytical results with experimental data obtained from tests performed at the Lewis 
Research Center with a J85-13 engine. In addition, comparison was also made with 
performance data obtained from the engine manufacturer. The comparisons established 
the validity of the simulation technique. 


INTRODUCTION 


The requirements in supersonic aircraft for high propulsion system performance 
can result in normal engine operation close to compressor stall and inlet operation 
close to unstart. To assure satisfactory overall operation under such conditions, pro- 
pulsion system design not only must include the dynamics of the inlet, engine, and asso- 
ciated controls but must also include the dynamic interaction between these system 
components. Consideration should be given to propulsion system design problems in 
the preliminary phases of any development program. Mathematical modeling and sys- 
tem simulation are generally used to resolve propulsion system design problems 
throughout a development program. In this way, expensive wind-tunnel and flight- test 
programs can be minimized, and system problems may be avoided or minimized. This 



report presents the development of mathematical modeling and simulation techniques 
for jet engine representation in an overall propulsion system study. 

Turbojet engine simulations have been developed by both eqgine manufacturers and 
government agencies. Generally, however, the models use linear transfer function 
techniques for system dynamics and overall performance maps for component steady- 
state behavior. Such representations of system dynamics are satisfactory for engine 
controls development and low-frequency (10 Hz or less) system analysis. But inlet 
flow and shock position dynamics tend to contain significant dynamics at higher frequen- 
cies. Engine representations for use in overall propulsion system studies should thus 
be valid for frequencies up to 50 hertz or more. The purpose of the study was to apply 
recently developed (refs. 1 and 2) high-frequency compressor modeling techniques to a 
complete jet engine model. In this way a simulation model more suitable for overall 
propulsion- system analysis than the linear transfer function approach could be estab- 
lished. 

This report proposes and develops simulation techniques applicable to any air- 
breathing engine. The jet engine simulation is capable of satisfying the following re- 
quirements: (1) representing the engine operating range from idle to rated speed, 

(2) representing engine dynamics to at least 50 hertz, and (3) representing engine opera- 
tion accurately at operating conditions off the normal engine operating line. The 
methods developed are applied to a J85-13 engine. Results obtained from the simulation 
are compared with experimental data. Agreement is demonstrated between the simula- 
tion and the experimental data to frequencies as high as 60 hertz. 


ENGINE DESCRIPTION 

The simulation techniques developed in this report are applied to the J85-13 engine. 
The J85 is an afterburning, turbojet engine with an eight- stage axial- flow compressor 
and a two- stage turbine. Compressor discharge air that bypasses the combustion zone 
is used to cool combustor and turbine components. The engine, for which a cutaway 
section is shown in figure 1 also has variable inlet guide vanes, controlled interstage 
bleeds, and a variable- area exhaust nozzle. Engine control is provided through hydro- 
mechanical main and afterburner fuel controls and the variable-area exhaust nozzle. 

The engine simulation included the compressor variable geometry feature s"but did 
not include afterburner combustion. The gas dynamics for the afterburner volume, 
however, were included. 
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ANALYTICAL MODEL 


For a dynamic simulation to be valid over a wide operating range, the simulation 
must exhibit realistic steady- state behavior over the entire operating range. Compo- 
nent steady- state operation can be determined from steady- state maps expressed in 
terms of various flow and thermodynamic parameters. The data for such maps are, in 
general, obtained from experimental data, although, in the early phases of a program, 
it may have been obtained from design studies. 

In this study system flow dynamics are included with the steady- state component 
maps without changing the steady- state operating conditions. The dynamics are included 
in such a manner as to separate the component maps and help avoid algebraic loops or 
iterations that might otherwise be encountered. Where possible, the system dynamics 
are formulated in terms of the same variables used in the component maps, thus simpli- 
fying the overall simulation. 


GENERAL GAS DYNAMIC EQUATIONS 

The basic representations used for the system gas dynamics can be derived from 
conservation equations. For quasi- one- dimensional flow the basic continuity, momen- 
tum, and energy equations can be written as 


— (p A) + — (p Av) = 0 
9t s 3x s 

(i) 

— (p Av) + — (p Av 2 ) - - Ag 0Ps 

at s 3x s ax 

(2) 

j-(p s Au t )+^- (p AvH) = 0 

at ax 

(3) 


The symbols used in these equations and throughout the text are defined in appendix A. 
The conservation relations in the form of equations (1) to (3) do not include friction 
losses, work done on or by the fluid, or the effects of external heat. However, these 
effects are in general included in the component maps. 

In the transition from the partial- derivative distributed- parameter form (eqs. (1) 
to (3)) to the total- derivative lumped- parameter form used in system simulation, the 
working variables must be chosen. Compressor and turbine maps are generally based 
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on total pressures and temperatures and on flow; thus, these variables were selected 
as the basic system parameters. Although some approximations (discussed later iri the 
text) were used in the translation from static to total quantities, in general, this choice 
of system variables was satisfactory. 

It is shown in appendix B that equations (1) to (3) can be approximated in lumped 
parameter form by 


- Oy = - - aw (4) 

dt s V 


A (w) = - AP 

dt l 


A (p T,) = - 1- A(T t W) 

m & L tr 


(5) 

(6) 


Equations (4) to (6) and an equation of state form the basis for the lumped- volume gas 
dynamics used in the simulation. 


COMPRESSOR 

Description 

The J85 compressor consists of eight stages with variable inlet guide vanes and 
interstage bleed flow. The first- stage pressure rise and flow characteristics and, in 
turn, the overall compressor map are thus functions of inlet guide vane angle. The 
overall compressor map is also a function of bleed flow. Such variable geometry fea- 
tures can greatly complicate compressor simulation. The stage stacking approach, 
however, can be used to avoid some of this complexity. 


Schematic Model 

Compressor modeling can be achieved by several known techniques. A proven 
method involves the simulation of an overall compressor map. Such maps can provide 
pressure, temperature, and stall information as a function of inlet parameters and en- 
gine speed. However, this technique provides only steady- state performance and, thus, 
does not consider interstage gas dynamics. 
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Figure 2 . - Schematic of n* h compressor stage. 

Another method of representing multistage axial- flow compressors is based on in- 
dividual stage modeling and on computing the pressure and temperature rise across each 
stage. These stages are interconnected, or stacked, to form the compressor model 
where the discharge conditions of one stage are used as the inlet conditions to the follow- 
ing stage. The gas dynamics are readily included by considering a volume associated 
with each stage. 

The stage stacking technique has several advantages over other simulation methods. 
The main advantage is that only functions of one variable must be considered. Each 
stage can be represented by a pair of pressure and temperature rise maps followed by 
its appropriate stage volume. Continuity, momentum, and energy equations are then 
used to compute the dynamic characteristics of each stage. A second advantage is that 
variable inlet geometry and interstage bleed can easily be included in the simulation. 

In addition, with the inclusion of the gas dynamics, the compressor simulation becomes 
unstable for certain operating conditions. The locus of these critical pressure ratios 
matches the actual compressor stall line (ref. 2). The disadvantage of this simulation 
technique is that a relatively large analog computer facility is required. Also, because 
of the large amount of nonlinear equipment, such as multipliers and function generators, 
the overall steady- state accuracy of the compressor model is limited. 

An idealized stage is shown schematically in figure 2. In essence, the basic 
assumption is that each compressor stage consists of a stage and its associated volume. 
The mathematical analyses for the static and dynamic performance are presented in 
appendixes C and D. 


Mathematical Model 

The stage stacking technique just discussed was used for the J85 compressor simu- 
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lation. As mentioned previously, this method is suited for this application because it 
readily permits the inclusion of the gas dynamics to study compressor stability. The 
compressor stage characteristics are presented as nondimensionalized pressure and 
temperature coefficients as functions of a flow coefficient. These coefficients are de- 
fined by the relations 


w-P - C n T tv,n-l 
“ N 2 


tc, n 
L\ P tv, n- b 


v2/7 


- 1 


^n 2 ( T tc,n “ T tv,n-1^ 
NP 


<P 


*n 


c,n jq- zc,n 


(7) 

(8) 

(9) 


where 



( 10 ) 


K n = 


a 


nr 


n 


(ID 


The axial component of the stage inlet air flow v z appearing in equation (9) must be 
determined to compute the flow function cp . This relation is derived in appendix C as 


v 

zc,n 


W c,n RT sv, n-1 

A P 
c,n sv, n-1 


(12) 


The pressure and temperature for equation (12) are defined in static parameters. Be- 
cause the compressor simulation is in terms of total parameters, an equivalent expres- 
sion for v„ was derived in appendix C as 
z 

^c,nV^v,n-l _ v zc,n 
5 v, n- l A c , n |/ 0 v,n-l 



7 



The dynamic characteristics of the compressor were simulated by assuming an in- 
dividual volume for each compressor stage. This volume was considered to be down- 
stream of each stage. The dynamic characteristics were obtained by applying continuity, 
momentum, and energy equations. The relations used for the gas dynamics were derived 
in appendix D as 


d Vn>=^c,n-W C; n + l-\ n > 


dt 


n 


L n 


dt ^sv, n T tv, n^ " v < T tc, n W c, n " T tv, n W c, n+1 “ T tv, n W b, n^ 
v n 


11. 

The equation of state for the n volume is as follows : 


(14) 


(15) 

(16) 


P tv, n R ^sv, n T tv, n 


(17) 


The interstage bleed was computed for the third, fourth, and fifth stages by assuming 
that the bleed ports were choked. The bleed flow can be computed from the relation 


W 


= k l A l 


tv, n 


b,n b b,n /= 

V T tv, n 


(18) 


Simulation Model 

The eight stages of the J85 compressor were simulated by using equations (7) to (11) 

and (13) to (18) combined with the appropriate pressure and temperature maps. The 

compressor stage maps are presented in terms of a pressure coefficient and a 

T 

temperature coefficient ip as functions of a flow coefficient cp . Typical curves for 
the second stage of the J85 compressor are presented in figure 3. These curves must 
be entered through the flow function so that a single-valued output results. The flow 
function cp can be computed from equations (9) and (13). 

The air angle 0, which appears in equation (13), varies as a function of corrected 
speed. The angles were assumed to be constant over the normal operating range of the 
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Figure 3. - Pressure and temperature maps for 
second stage of compressor. 


J85 engine because the deviations from average values were small. 

Equation (13) was solved on a digital computer for v z as a function of inlet air flow 
and angle ft. The results were approximated with a quadratic expression having a stand- 
ard error of estimate of 2. 93 meters per second (9. 62 ft/sec). A typical curve for the 
analog simulation is presented in figure 4. 

The stage maps for the first compressor stage must be modified to account for the 
effects of inlet guide vane variations. First- stage pressure and temperature character- 
istics were generated for closed, half open, and full open inlet guide vane positions, for 
which the pressure maps are shown in figure 5. The pressure and temperature charac- 
teristics were then determined by interpolating between these curves. A detailed dis- 
cussion of this and several other techniques for bivariant function generation is pre- 
sented in appendix E. 

Interstage bleed was included in the third, fourth, and fifth compressor stages. "-The 
choked flow equation (18) was used to compute the actual bleed flow. The area of the 
bleed ports was scheduled as a function of corrected speed and compressor inlet temper- 
ature. The bleed valve movement occurs within the 80 to 100 percent of the corrected 
speed range. 
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Air velocity, v z /\/ 6, ft/ sec 



Corrected weight flow, W VG/ (6A), lbm/(ft 2 )(sec) 
Figure 4. - Typical air velocity computation. 


Inlet guide vane position 



Normalized flow coefficient, (pl<p 


Figure 5. - First stage pressure map for various inlet guide vane positions. 
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The outlet guide vanes were represented by a volume and pressure drop. The tem- 
perature within the guide vane volume was assumed to be that of the compressor dis 
charge. 

A block diagram of a typical compressor stage showing the computational techniques 
used in the simulation is presented in figure 6. As observed, the stage inlet conditions 
are used to compute the value of the flow function. This value was then used to deter- 
mine the stage pressure and temperature coefficients from which the discharge condi- 



\ n+l 


Figure 6. - Block diagram for simulation of compressor stage dynamics. 


tions were computed. The stage discharge conditions combined with the gas dynamics 
are used to compute the pressures and temperatures of each stage volume. These 
parameters were then used for the inlet conditions to the next stage. 

The pressures and temperatures in the compressor simulations are computed in 
terms of total parameters. The simulation neglects the difference between static and 
total parameters because of the complexity involved in computing the Mach number for 
each stage. Further discussion is presented in appendix D and reference 2. 
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COMBUSTOR 


Description 

The J85 annular combustor is shown in figure 7. Air from the compressor enters 
the liner through perforations (thimble holes) and louvers. The thimble holes direct air 
into the burning area (primary combustion flow), while the louvers provide a boundary 

layer (secondary cooling flow) of comparatively cool air along the inner surface of the 
liner. 



Figure 7. - J85-13 combustor. 
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Performance 


In modeling the combustor steady- state operation pressure drop, heat addition and 
combustor efficiency must be accounted for. 

Pressure drop . - The total pressure drop across a combustor can be attributed to 
fluid friction and momentum loss due to heat addition and stream mixing. The analytical 
representation is based on the assumption that the combustion chamber can be divided 
into equivalent fluid friction and heat addition zones. The air enters the combustion zone 
and experiences a total pressure drop due to friction. This pressure drop corresponds 
to the actual mixing drop that occurs without combustion, which can be determined ex- 
perimentally and correlated with theory. Momentum pressure loss is assumed to occur 
in the combustion zones with the addition of heat to the air. This drop is assumed to 
occur after the friction loss. The overall loss in pressure is then formed from the sum 
of the two losses. The resulting expression is 

Wo 

AP t,4~( K A T t,3 + K B T t, 4 ) (19) 

\ 3 

where K^ and Kg are experimentally determined coefficients peculiar to each com- 
bustor: is found from noncombustion flow tests, and Kg from combustion flow 

tests. 

Heat addition . - The actual combustion involves processes such as atomization, 
vaporization, mixing, and chemical reaction. These processes are gradual in time and 
space. For combustor simulation, however, it is only necessary to determine relations 
between initial and final conditions. 

Since variations in specific heat due to composition and temperature changes can 
introduce significant errors, variable specific heat is considered in the combustor per- 
formance. However, it is assumed that 

h 4 = f < T t, 4> (20) 

Using data obtained from reference 3, this functional relation was evaluated and is ^ 
plotted in figure 8. 

Efficiency . - Because of the complexity of the combustion process, no exact analyt- 
ical treatment of combustor efficiency was attempted. For the purposes of this simula- 
tion, combustor efficiency was approximated by a correlation of the form 

( 21 ) 


^4 “ f ( P t, 4 A ^t, 4^ 
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Combustion enthalpy, h 4 , Btu/lbm 
Figure 8. - Temperature of combustion gas as function of enthalpy. 


The efficiency 77 4 is here defined as that portion of the heat of combustion that is avail- 
able for a gas temperature rise. The parameter P t 4 AT t 4 is based on ideal combus- 
tion temperature rise and includes the effects of fuel- to- air ratio and inlet pressure 
variations but neglects the effect of inlet temperature variation. 

Combustor efficiency for the J85 is shown in figure 9. Steady- state data on and 
near the operating line at sea level and altitude were used to generate the relation. It 
provides reasonable accuracy over a wide range of operating conditions. 

Schematic Model 

The combustor model used in this simulation is formed by adding dynamics*to the 
steady- state characteristics presented in the performance section. The combustor is 
lumped into a single equivalent one- dimensional volume. The volume is bounded by the 
inner and outer shells, with its length being defined by the compressor exit face and the 
turbine first- stage stator face. 
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Compressor Combustor 

Figure 10. - Combustor model. 


Turbine 


In the formation of the combustor model through the addition of dynamics to the ^ 
steady- state characteristics, provision is made for heat addition. The combustion pro- 
cess is assumed to be continuous in time and external to the basic combustor gas 
dynamic lumped volume. The end product of the combustion process is assumed to be 
heat, which is added to the lumped volume. Schematically, this representation is shown 
in figure 10. 
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Mathematical Model 


Applying the basic equations of momentum and energy to the schematic model re- 
sults in the following mathematical representation: for combustor weight flow the 
momentum equation yields 


where 


dW 4 gA 4 




>- AP t, 4 ] 


AP, 


t,4 


K C^1 

P t,3 


(0.771 T t s - 0.085 T t 4 ) 


( 22 ) 


(23) 


The constants in this equation are empirically determined from J85 turbojet experimen- 
tal data. 

The energy equation gives 


_d 

dt 


^s, 4*t, 4^ 


_V_ 

V, 


[W 4 T t,3 


< T t ,4 


W 5 T t,4 




(24) 


Combustion limits were not included in this representation. 

The fuel flow W f is the result of combustion dynamics imposed on the fuel flow 
external to the lumped volume gas dynamics. 


• t t 



Ke~ gs 

(1 + t a s )(1 + TgS) 


(25) 


where and r B are the time constants associated with heat release due to combus- 
tion lag. The delay time a is associated with the ignition lag and transport delay due 
to the difference between where the fuel flow is measured and where the fuel nozzle is. 
The continuity equation gives 


-^s,4> =^(W 4 + W; T - W 5 ) 
V 4 


dt 


( 26 ) 
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Simulation Model 


The mathematical representation was simulated on the analog computer in a straight 
forward manner. 

The burner efficiency curve shown in figure 9 was approximated using a diode func- 
tion generator. The fuel flow dead time was generated using a fourth- order Pad6 
approximation. 


TURBINE 

Description 

The J85-13 engine uses a two- stage turbine cooled by secondary combustion air. 

The most direct approach to turbine modeling would have been to apply a stage stacking 
technique, similar to that used for the compressor, to the two turbine stages. Individual 
stage data for the turbine were not available, however, and it was thus necessary to 
represent the turbine from overall performance data. 


Schematic Model 


The turbine data available during the J85 simulation program was in the form of an 
overall performance map. These data are plotted in figure 11. The direct application 
of information in this form would be extremely difficult on an analog computer. To 
avoid this problem the data were manipulated to a form more amenable to analog simu- 
lati on t e c hnique s . 

The data in the performance map of figure 11 can be regarded as giving a turbine 
enthalpy parameter and turbine efficiency as functions of flow and speed parameters; 
that is, 


Ah 5 * /V T t, 4 

nn = *(l)l p ’ 
T t, 4 \ *t,4 



^5 “ f (2) 


rW 4l^4 


t, 4 



(27) 


(28) 
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Figure 11. - Overall turbine performance map. 
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The turbine efficiency tj 5 can be defined as 


^5 " 


Ahg 

Ah' 5 


(29) 


where Ah£ is the ideal isentropic enthalpy drop. Using the isentropic relation between 
pressure and temperature, the pressure ratio across the turbine can be expressed as 




Ah 5 y/^-o 
v 5 c p T t, J 


Considering equations (27) to (30), it can be concluded that 



(30) 


( 31 ) 
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These functions are shown in figures 11 and 12. The enthalpy parameter is superim- 
posed on the overall performance map. 


Turbine pressure ratio, 
P t,5 /P t,4 



100 150 200 250 300 350 400 450 

Speed parameter, rpm/°R 1/2 

Figure 12. - Turbine flow performance. 


The functional relations of equations (31) and (32) were normalized with a turbine 
speed parameter to separate and make them well behaved for the simulation. The re- 
sulting curves are plotted in figures 13 and 14. In this representation, the enthalpy func- 
tion becomes a torque parameter, that is, 
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Figure 13. - Turbine torque characteristics. 









and the normalized flow parameter becomes 


W 5 T t,4 ,„ W 5_ 


NP t, 4 Np s, 4 



( 34 ) 


An attempt at a correlation leading to a function of a single variable is presented in 
appendix F. The results of this correlation were good, but not accurate enough for this 
case to duplicate known steady- state data. Therefore, the two- variable correlation was 
used. 

No attempt is made to model the actual turbine cooling flow process. The flow is 
instead added to the appropriate volume and mixed. The momentum of the stream is 
neglected. Since performance is based on inlet flow, no appreciable error is introduced. 
The turbine model is shown schematically in figure 15. 


Combustor 


Turbine maps and steady- 
state calculations 


Turbine volume 
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Z T ; 



%; T t, 3 

Figure 15. - Turbine model. 


Mathematical Model 


The turbine is represented mathematically by the imposition of the basic conserva- 
tion equations on the turbine steady- state performance characteristics. Input variables 
and concurrent output values are used to calculate steady- state values of the enthalpy 
and flow parameters. The inertance is extremely low for this model, thus, the turbine 
flow was determined directly from the steady- state performance map 


W 5 T t,4 _ f 


NP. 


t, 4 



N 


,4 V T t, 


(35) 


given in figure 14. 
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The turbine enthalpy drop is computed from 


Ah 5 




(36) 


which is plotted in figure 13. The enthalpy drop determined in equation (36) is used to 
compute the turbine discharge temperature from 



(37) 


The pressure and temperature downstream of the turbine are determined from con- 
servation equations written for the volume. The continuity equation yields 

^<M=^(W 5 + W tc -W 8 ) (38) 

8 

The energy equation is 


i (p s, 8 T t, 8> =f g < T t, 5^5 - T t> 3 W tc - T t> 8 W 8 ) 


(39) 


AFTERBURNER 

The afterburner volume was included in the simulation; however, afterburner com 
bustion was not considered. The lumped volume representation consisted of mass and 


Turt) ine Afterburner 

volume volume 

dynamics dynamics 
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p u 


*5.1 


OO 

T t, 8 ( 





1 » 


Nozzle 



Figure 16. - Lumped volume dynamics. 
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momentum balances, but no energy balance. The lumped volume dynamics are shown 
in figure 16. 

The momentum and continuity equations are 



) = 


-5.1 g 

*5.1 




(40; 


jl 

dt 


( P t,8)=^-( W 5.1 T t,8- W 8 T t,8) 

v 8 


(41) 


Inherent in this representation is the assumption that the afterburner volume is at the 
same temperature as the turbine discharge. 


EXHAUST NOZZLE 
Schematic Model 

The nozzle was simply considered as a variable-area flow passage capable of chok- 
ing. It does not, as such, have an associated volume and becomes essentially the termi- 
nal boundary of the turbojet. 


Mathematical Model 


The nozzle is represented mathematically by the relation 


W ( 




t, 8 



(42) 


where P g Q is ambient static pressure. The parameter K g is variable and lumps the 
nozzle flow coefficient and nozzle area with other fluid constants. Thus, it is propoi> 
tional to the nozzle area and represents the turbojet engine terminal impedance. 

The function f(P g Q /P t g ) is the well-known compressible flow function: 


f 





(43) 
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For pressure ratio values below critical, it becomes constant and simulates a choked 
passage. 


Simulation Model 

The compressible flow function was approximated by a fixed diode function gener- 
ator. The parameter Kg was scheduled as a function of engine speed to match the 
steady- state operating line. 

ROTOR DYNAMICS 

The operating line of a turbojet engine is the locus of steady- state performance 
match points between the compressor and turbine. A mismatch between these compo- 
nents produces an unbalanced torque or acceleration potential, which must be integrated 
through a dynamic relation to change the speed and seek a steady- state match. 

The relation governing rotor dynamics is the conservation of angular momentum: 

X =E L ( 4 4) 


Using the expression 


(P = o)L = JW Ah 

the dynamic relation between engine speed and enthalpy change is 


(45) 


dN = /30\ 2 X E WAh 

dt W 


(46) 


The net sum of enthalpy differences are those that occur across the compressor and 
turbine, thus — 


dN 

dt 



5 h 4 " ^8 h 5^ ~ ^3 h 3 + ^b h b ~ 


w 2 ty] 


(47) 


where h fo is compressor bleed total enthalpy. 
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SIMULATION VERIFICATION 


To establish the validity of the engine simulation described in this report, data ob- 
tained from three sources were compared with results obtained from the simulation. 
First, steady- state compressor operation and the engine system normal operating line 
were compared with data obtained from the manufacturer’s digital engine program. 
Second, compressor steady- state and stall line data obtained from reference 4 were 
used to check the compressor model. Finally, unpublished engine system frequency 
response data were used to verify the engine system simulation’s dynamic response. 

Data for compressor operation along constant speed lines and the compressor stall 
line were obtained from a test program run at Lewis. The test results were obtained for 
an off- schedule inlet guide vane and inlet bleed configuration. Details of this test pro- 
gram can be found in reference 4. 

Engine system frequency response data were obtained from a test program con- 
ducted in the 10- by 10- foot Supersonic Wind Tunnel. The program used a J85-13 en- 
gine and a mixed- compression supersonic inlet. Details of the inlet dynamics are pre- 
sented in reference 5. The engine system was instrumented to obtain total and static 
pressures at various engine and inlet stations. Frequency response data were obtained 
by varying combustor fuel flow and recording system pressures. The frequency range 
of available data is limited to about 60 hertz. 


STEADY-STATE PERFORMANCE 
Compressor Mapping 

To separate the compressor from the combustor and turbine for independent analy- 
sis, the compressor model was terminated with an equivalent volume and a choked ori- 
fice. This simplification was used during the mapping of the steady- state performance 
characteristics of the compressor. 

The compressor performance map can be generated by varying the terminating im- 
pedance at various rotational speeds. With increasing compressor terminating imped- 
ance, caused by reductions in the choked station area, the compressor operating point 
moves along the constant speed line. Several constant speed lines were mapped. A- 
representative compressor map, for a limited speed range, obtained from the simulation 
is presented in figure 17. As a comparison, sea- level experimental data obtained from 
reference 4 are also presented in this figure. Both the simulation and test results were 
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Corrected airflow, Ibm/sec 
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speed, 

N/yf, 
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Figure 17. - Compressor map for off -schedule inlet guide vanes and interstage bleeds. 


obtained for an off- schedule inlet guide vane and interstage bleed configuration. The 
agreement between the experimental and analytical results is good and verifies the 
steady- state representation of the compressor. 


Compressor Stall 

As described in reference 2, it was found that, as the terminating impedance was 
increased, the simulation will exhibit a nonlinear oscillation at an operating point on the 
constant speed line. Further increase in the terminating impedance causes a higher 
amplitude oscillation. The simulation exhibited hysteresis in that the impedance must 
be reduced to a level below the point of initial instability to reestablish system stability. 

The locus of the operating points just before sustained oscillations is defined as the 
compressor stall line and is shown in figure 18. Again as a comparison, the experimen- 
tal stall line is also presented. The figure indicates that the analytical and experimental 
stall lines compare favorably. 

Some understanding of compressor performance in the vicinity of stall can be ob- 
tained from the individual stage pressure maps. The J85 normalized stage characteris- 
tics for corrected speeds of 100 and 80 percent are illustrated in figures 19 and 20, re- 
spectively. Only the first- stage characteristics are varied with the inlet guide vane 
position. 
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O Normal operation 
A Near stall 



speed. 


The normal operating points of the various stages (equilibrium point) are indicated 
on the compressor stage maps. The effects of the variable geometry is evident in that 
the weight flow through the latter stages of the compressor is reduced at the lower 
engine speed. Interstage bleed increases the weight flow through the front compressor 
stages so as to avoid stall. With higher weight flow, the flow coefficient <p increases 
to shift the operating point to the stable region. As is also observed from these curves, 
at higher operating speeds most of the compressor pressure rise is achieved in the lat- 
ter stages. At lower operating speeds, the pressure rise is obtained in the front stages. 

As can be further observed, normal operation at the design speed occurs on the neg- 
ative slope portion of the pressure- rise maps for all stages. As the stability boundary 
is approached, all the stages move towards a less negative slope. The second and third 
stages move into the positive slope region. For operation at 80 percent of design speed, 
normal operation of stages 2 and 3 occurs within the positive slope region. As the sta- 
bility boundary is approached, these two stages shift towards increasing positive slope. 

After the stability boundary is crossed, the simulation is invalid because the stage 
characteristics are inaccurate for stalled operation. Any degradation of the stage per- 
formance in the stalled region is not considered in this engine simulation. A detailed 
discussion of compressor stall can be found in reference 6. 
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TABLE I. - STEADY- STATE PERFORMANCE 


Design engine 
corrected 
speed, 
percent 

Type of data 

48 

Experimental 

Analytical 

80 

Experimental 

Analytical 

100 

Experimental 

Analytical 


Compressor discharge 


Pressure Temperature 


Combustor 


Pressure Temperature 


kN/m 2 I lbf/ft 2 K °R kg/sec lbm/sec kN/m 2 Lbf/ft 2 


64. 8 3 442 341 


0 424 763 


9.6 14 820 550 990 20.0 


6.21 

■ 

6.08 

1 

.3.9 

3 

.4.5 

3 


53. 1 13 640 

53.1 13 640 


837 1507 
853 1535 


kg/sec lbm/sec kN/m lbf/ft' 


4.8 2189 

5.5 2203 


Turbine 


Pressure Temperature Air flow 


K °R kg/sec lbm/sec 


2206 19.2 

2232 19.2 


KSlISiS 


647 1164 
643 1158 


42.4 246.2 5141 


992 1785 19.6 


42.4 264.1 5515 1021 1838 19.6 




























































































Engine Performance 


To verify the complete engine simulation over its entire operating range, analytical 
values of pressure, temperature, and weight flow obtained from the simulation were 
compared with experimental values. They were also compared with steady- state data 
obtained from the manufacturer’s overall digital performance program. The results are 
shown in table I. 


DYNAMIC PERFORMANCE 

This section presents a comparison of frequency response obtained from the engine 
simulation with experimental data obtained during actual engine tests in the Lewis 10- 
by 10- Foot Supersonic Wind Tunnel. During the tunnel test program an axisymmetric 
mixed- compression inlet was used with the engine. For comparison with the tunnel test 
data, a linear inlet representation based on the techniques of reference 7 was added to 
the engine model. The inlet representation is only valid for small perturbations so that 
the inlet was modeled for the designated engine operating conditions. The inlet was re- 
quired to assure correct compressor face conditions. 

In the actual test program, a long feed line was used between the high- response 
fuel- control valve and the engine fuel spray nozzle. Since this line introduced unwanted 
dynamics to the system, it was necessary to modify the experimental data and eliminate 
these dynamics. The frequency response curves presented in this section are therefore 
all referenced to the combustor fuel spray nozzle pressure. 


Rotor Speed Effects 

The computer representation of the J85 engine was developed to establish the capa- 
bility of modeling engine system responses in the frequency range from 10 to 50 hertz. 
Engine system low-frequency performance is dominated by the rotor- speed - time con- 
stant. Since this time constant has a break frequency below 1 hertz, it was assumed that 
unpublished test data could be matched in the frequency range of interest (above-10 Hz) 
by operating the engine simulation at constant speed. The frequency response data pre- 
sented in this report for the engine simulation were thus obtained by operating the com- 
puter model in a constant- speed mode. 

To establish the validity of the assumption that rotor dynamics could be neglected 
in the frequency range above 10 hertz, a brief linear analysis was performed. The re- 
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sponse of combustor total pressures to fuel flow disturbances both with and without rotor 
speed effects was determined and compared. Conceptually, the response of combustor 
total pressure to fuel flow disturbances can be approximated by 




3P t, 4 t AN 8P t,4 \ 
3wj' AWf 3N j 


(48) 


The first term of equation (48) represents the rise in combustor total pressure, which 
is caused by energy addition from combustion. This rise would include combustion dy- 
namics and is characteristically fast in response. The second term represents the in- 
crease in pressure due to engine rotor dynamics, which is dominated by the rotor time 
constant and is thus slow in response. An approximate transfer function for this overall 
process can be given as 


P t,4, G ;5A . N ^, ( K r + K R ) 


1 + 


( K G t R + k r t g^ s 

k g + k r 


w. 


w x 


(1 + r R s)(l + 7 a s) 


(49) 


G" 


where r R and r Q are the rotor and gas time constants, respectively. The combustor 
total- pressure frequency response, neglecting transport delays, can be determined from 
equation (49). The results of such a calculation using t r = 0. 5 second, Tq - 2.9 milli- 
seconds, and appropriate gains are shown in figure 21. Also plotted in figure 21 are the 
results obtained from equation (48) where the rotor dynamics are neglected. As can be 
seen, the effects of rotor dynamics on combustor total pressure response are significant 
only below 5 hertz. It should be noted, however, that in the range from 5 to about 15 
hertz some phase lag (7° to 2°) is added by the rotor dynamics. The use of constant 
rotor speed in the simulation is thus a reasonable assumption in the determination of 
system responses above 10 hertz. 


Combustor 

The frequency response of combustor total pressure to spray nozzle pressure per- 
turbations is presented in figure 22. Appropriate values for ignition delay and combus- 
tion lag were obtained from the experimental test results of reference 8 for inclusion in 
the simulation. As a result, the dynamic response of the combustor exit pressure of the 
simulation was closely matched to the engine experimental data. 
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Compressor 


The frequency response of compressor discharge total pressure to spray nozzle 
pressure is shown in figure 23. The attenuation characteristics are nearly identical; 
however, a small additional phase shift is exhibited by the actual J85 engine data. 

The frequency response of compressor discharge static pressure to spray nozzle 
pressure is presented in figure 24. The conversion from total to static pressure was 
accomplished with a linearized conversion formula based on pressure and flow that was 
obtained from reference 7. 

The frequency response of the sixth- stage static pressure to spray nozzle pressure 
is presented in figure 25. The normalized frequency response curves for the sixth stage 
are quite similar to the compressor discharge responses when they are perturbed with a 
fuel flow disturbance due to the small volume between these stages. 

The amplitude of upstream compressor stage pressures are less sensitive to fuel 
flow disturbances than the downstream stage pressures. This trend was exhibited by the 
third- stage response. The amplitude of the pressure variations become extremely small 
near the compressor face, thus inaccuracies in the test data can exist. The response of 
compressor third- stage static pressure to spray nozzle pressure is shown in figure 26. 

It should be noted that there is some deviation between the analytical and experimental 
results in the low frequency region. 

It should be noted that, in general, the phase shift of the actual J85 engine is slightly 
higher than that obtained for the engine model. The slight difference can be attributed 
to the rotor dynamics effect at low frequencies. The analytical phase shift of combustion 
to spray nozzle pressure (fig. 21(b)) indicates that additional phase shift is obtained from 
the rotor dynamics. It is probable that if the speed effect had been included during the 
frequency- response determinations, better agreement between experimental and analyt- 
ical results could have been obtained at frequencies below 10 hertz. 


SUMMARY OF RESULTS 

Overall, the results of this study indicate the validity of the simulation over the 
range investigated, considering, of course, the accuracy of the experimental data to 
which it was compared. The techniques used in the simulation can be applied to other 
engines to obtain reasonably accurate steady- state and transient representations. 

The study confirms that a mathematical representation using the dynamics inherent 
in the conservation equations and engine geometry will provide a better simulation than 
those using component representation and linearized dynamics. The flexibility and wide- 
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range nonlinear capability outweigh the large equipment complement. The dynamic stage 
stacking technique reproduces the experimental steady- state compressor operating map 
very well and provides a method of predicting the stall line as well. This technique is a 
strong analytical tool especially for engine studies where the stall line is not known. 

The technique also provides the capability to couple inlet dynamics to the engine and 
predict propulsion system interaction. 

The combustor representation seems to be the least rigorous in the engine simula- 
tion. The difficulty lies in obtaining a simple, but complete, representation of the com- 
plicated process. 

The results also indicate that care should be exercised in assuming constant speed 
for simulated dynamic response. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 24, 1971, 

132-15. 
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APPENDIX A 



SYMBOLS 


A 

2 n.2 

area, m ; ft 

k d 

coefficient, (kg)(m 2 )(K 1//2 )/ (N)(sec); 

a 

speed of sound, m/sec; 

(lbm ) (ft 2 ) (°R 1/2 )/(lbf ) (sec) 


ft/sec 

k g 

coefficient, (N)(sec)/(kg)(m 2 ); (lbf) 

C n 

coefficient, rpm 2 /K; 

(sec)/(lbm)(ft 2 ) 

rpm 2 /°R 

K n 

coefficient, (rpm)(sec)/m; (rpm) 

C P 

specific heat at constant 


(sec)/ft 

pressure, J/(g)(K); 

k r 

coefficient, (N)(sec)/(kg)(m 2 ); 

c v 

Btu/(lbm)(°R) 
specific heat at constant 

K 8 

(lbf ) (se c )/ (lbm ) (f t 2 ) 
nozzle parameter, (kg)(m 2 )(K 1//2 )/ 

g 

volume, J/(g)(K); Btu/ 
(lbm)(°R) 

gravitational constant, 

k b 

(N)(sec); (lbm)(ft 2 )(°R 1//2 )/(lbf) 
(sec) 

bleed flow coefficient, (kgXK*/ 2 )/ 
(N)(sec); (lbm)(°R 1/,2 )/(lbf)(sec) 

torque, N-m; ft- lbf 


l(kg)(m)/(N)(sec 2 ); 32.17 
(lbm ) (ft )/ (lbf ) (s e c 2 ) 

L 

H 

total enthalpy, J/g; Btu/lbm 

l 

length, m; ft 

h 

static enthalpy, J/g; 

M 

Mach number 

Ah 

Ah’ 

Btu/lbm 

actual isentropic expansion 
value, J/g; Btu/lbm 
ideal isentropic expansion 

N 

N f 

n 

rotational speed, rpm 
nonlinear rotational speed, rpm 
number of mass particles 


value, J/g; Btu/lbm 

P 

pressure, N/m , lbf/ft z 

h c 

heat of combustion, J/g; 


power, N- m/sec; ft-lbf/sec 


Btu/lbm 

R 

universal gas constant, 287(N)(m)/ 

I 

moment of inertia, N-m- 


(kg)(K); 53.3 (lbf)(ft)/(lbm)(°R) 


sec 2 ; lbf- ft- sec 2 

r 

mean radius, m; ft 

J 

mechanical equivalent of 

r rp 

tip radius, m; ft 


heat, 1 N-m/J; 778.3 
ft- lbf /Btu 

X 

s 

Laplace operator, sec" * 

K 

coefficient 

T 

temperature, K; °R 

K A’ K B> k c 

coefficients, (N 2 )(sec 2 )/ 

AT 1 

ideal total temperature drop across 

(kg 2 )(m 4 )(K); (lbf 2 )(sec 2 )/ 
(lbm 2 )(ft 4 )(°R) 

t 

turbine, K; R 
time, sec 


40 



U mean rotor speed, m/sec; ft/sec 

Urp tip rotor speed, m/sec; ft/sec 

u internal energy, J/g; Btu/lbm 

V volume, m^; ft** 

v velocity, m/sec; ft/sec 

v^ tangential velocity, m/sec; 

ft/sec 

v_ axial velocity, m/sec; ft/sec 

z 

v zc axial velocity associated with 

stage characteristics, m/sec; 
ft/sec 

W weight, kg; lbm 

W weight flow, kg/sec; lbm/sec 

• 1 1 

W f externally acted upon fuel flow, 

kg/sec; lbm/sec 

x distance, m; ft 

y distance, m; ft 

z distance, m; ft 

z* distance, m; ft 

a coefficient, 30 sec/min 

j. 3 rotor air inlet angle, deg 

y ratio of specific heats, 1. 4 

6 ratio of total pressure to stand- 

ard atmospheric pressure 

?] efficiency 

9 ratio of total temperature to 

standard atmospheric temper- 
ature 

X work- speed parameter 

O O 

p weight density, kg/m ; lbm /ft 

or delay time, sec 

t a , r-g time constant, sec 


cp flow coefficient 

p 

pressure coefficient 
T 

temperature coefficient 
a ) angular velocity, sec * 

Subscripts: 

b variable associated with stage 

bleed 

c variable associated with stage 

characteristics 

f fuel 

G gas 

(i) denotes function, i = 1, 2, 3, ... 
in inlet 

N speed 

n stage number designation 

out outlet 

p variable associated with particle 

R rotor 

s static condition 

sc static condition variable associ- 

ated with stage characteristics 

stg 1 first stage 

stg 2 second stage 

sv static condition variable associ- 

ated with stage volume 

TC turbine cooling ^ 

t total condition 

tc total condition variable associated 

with stage characteristics 
tr total condition reference state 

tv total condition variable associated 

with stage volume 
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v variable associated with stage vol- 

ume 

0 free- stream condition 

2 compressor inlet 

3 compressor discharge 


4 combustor 

5 turbine 

5. 1 turbine discharge 
8 nozzle 
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APPENDIX B 


BASIC GAS DYNAMICS REPRESENTATION 

The relations used for the representation of system gas dynamics can be derived 
from the continuity, momentum, and energy equations. The mathematical manipulations 
required to formulate these equations in terms of the chosen variables are given in this 
appendix. 

The laws of conservation of mass, momentum, and energy can be written for quasi- 
one- dimensional flow as 


— (p„A) + JL (p Av) = 0 (Bl) 

3t S 0x S 


.1 (p Av) + A (P s Av 2 ) = - Ag — 5 (B2) 

at s ax s ax 


— (p Ail) + — (p AvH) = 0 (B3) 

at s 1 dx s 


Continuity equation: 

Noting that W = p c Av, equation (Bl) becomes 

b 


- (P g A) =-iw 

at s ax 


or 


2(p )..iiw 

at b a ax 


and, lumping in the spatial variable 


— (p s ) =- (W, - w 2 ) 
dt s V 


(B4) 


(B5) 


(B6) 
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Momentum equation: 

Substituting W = p 0 Av into equation (B2) yields 


a . a « 9P_, 

— (W) + -2- (p Av 2 ) = - Ag — ® 
St Sx s 3x 


O 

Neglecting the convective momentum term 9/9x (p_Av ) gives 

b 


a • SP„ 

— (W) = - Ag — § 
9t 9x 


or in lumped form 


_d 

dt 


( W )=^( P s,l- p S) 2) 


Energy equation: 

Substituting the term W = p g Av directly into equation (B3) gives 


— (P-Avl) + — (WH) = 0 
St s T Sx 


which becomes 


- Wt> - - “ 


at 


A 9x 


or 


_d 

dt 




t,i 


W 2 T t,2> 


(B7) 


(B8) 


(B9) 


(BIO) 


(Bll) 


(B12) 
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APPENDIX C 


STEADY-STATE COMPRESSOR STAGE REPRESENTATION 


The mathematical description of the steady- state compressor characteristics is 
based on the stage- stacking technique. This method uses two-dimensional maps to 
represent each stage. The individual stage performance is described by a pressure 
coefficient i j? and a temperature coefficient as functions of a flow coefficient <p . 
The coefficients are mathematically described by 

= (Cl) 

tt2 


* 


T 



(v z /Ve) 
(u /yfe) 


(C2) 


(C3) 


11 * 

The total enthalpy rise for the n stage can be related to the stage total pressure ratio 
from 


Ah tc, n c p T tv, n- 1 


p \(y-i)/r 
tc,n \ _ i 


tv, n- ly 


(C4) 


The adiabatic enthalpy rise can be computed from 

Ah tc, n " c p (T tc, n “ T tv, n- 1^ 
and the stage wheel speed U T from the relation 

?r Nr rp 

U T = i (C6) 

a 

The magnitudes of the inlet pressures and temperatures of the n 1 stage are determined 
from the exit conditions of the preceding stage. Therefore, equations (Cl) to (C3) can 
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be combined with equations (C4) to (C6) to develop a set of equations for the n th stage. 


* 


P _ C n T tv, n- 1 

n 2 
1ST 


p \2/7 
tc,n \ _ x 


tv, n- 1/ 


(C7) 


K ^( T tc,n- T tv,n-l) 
N 25 


(C8) 


_ K n 

^n N v zc,n 


(C9) 


where 


2gJc (a y 
C = 2 


n 


2 2 

* r T,n 


(CIO) 


K = 
n 


a 


irr 


(Cll) 


n 


The axial velocity v z can be computed from the flow, total pressure, and total 
temperature at the stage inlet. The flow and axial velocity are related by 


v z = 


W 

Ps A 


(C12) 


Eliminating the density from equation (C12) with the state equation, 


P s Ps RT s 


(C13) 


results in 


v z = 


WR T s 
A 


(C14) 
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or 


y =WR^(i + r^i 

z A P, \ 2 




M‘ 


Equation (C15) can be written in terms of corrected parameters as 






or 


A6 \ 2 


-l/(y-l) 


It can be shown from the energy equation that 


V 2 2Jc p L _ V 


2 yR 


T 4 


Therefore, 


T S_ 

1- y- 1 

/v\ 2 ' 

II 

i + y - 1 m 2 " 1 

T t 

CM 

1 

\V_ 


— i 
CM 


-1 


Combining equations (C17) and (C19) yields 


W >fe _ v z 
A6 


1 - 


v 


yfej 2gJc p T tr J 


l/(y-l) 


p tv 


The velocity term appearing inside the brackets of equation (C20) is the magnitude 
velocity vector. If the airflow approaches the rotor inlet at an angle 0, then 


v — v cos /3 


(C15) 

(C16) 

(C17) 

(C18) 

(C19) 

(C20) 
of the 
(C21) 
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and equation (C20) becomes 


wy'o = 


|l/(yi) 


l - 


W 2gJc p T tr cos 2 0 


aA ^0 

Equation (C22) can be rewritten for the n th stage as 


*tr 


W. 


:,n|/^v,n-l v zc, n 


A c, n^v, n- 1 n-1 


|i/(r-i) 


l - 


zc, n 


VVn-1/ 2gJc p T tr cos 2 ^ 


p tr 


(C22) 


(C23) 
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APPENDIX D 


STAGE GAS DYNAMICS 


The gas dynamics of each compressor stage are assumed to occur in an equivalent 
stage volume. This stage volume is located downstream of the stage pressure and tem- 
perature rise characteristics. The stage gas dynamics are represented by continuity, 
energy, and momentum balances across the volume. 

Assuming that no work is performed in the volume and that the net heat input is 
negligible, an energy balance can be written for the n^* 1 volume as 


V n dt ^sv, n u tv, n^ ^tc, n W c, n 



(Dl) 


Since u^. = c y T^. and h^. = c^T^., equation (Dl) becomes 


_d 

dt 


^sv, n T tv, n^ y 


< T tc, n W c, n " T tv, n W c, n+1^ 


n 


(D2) 


To determine the density of the gas in the stage volume, it is assumed that the rate 
of change of the static density is proportional to the net weight flow. 


_d 

dt 


^sv, n^ 



< W c,n- W c,n + 1> 


(D3) 


The total pressure in the volume can be computed from the equation of state: 


/ 2 \ 5/2 

P tv, n “ R r + 2M v, n) p sv, n T tv, n 


(D4) 


A momentum balance across the stage and associated volumes is used to determine 
the weight flow into each stage. The flow acceleration is caused by the difference be- 
tween the static compressor pressure rise AP 0 „ and the pressure rise 

11 

(P oir - P 1 ) imposed by the two stage volumes. This yields 
s v ; n sv ? n- i. 


_d 

dt 


A r,g 

(W c,n>=f- 

n 



sc, n 


(P. 


sv, n 


- P 


sv, n- 1! 


(D5) 
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However, the inlet conditions to the stage are equal to the outlet conditions of the pre- 
ceding stage volume, thus AP gc n = P gc n - P n-1 and equation (D5) reduces to 


A ( w 

dt 


A g 

=— (p c 


- P, 


(D6) 


If the Mach number at the stage exit and in the volume are identical, then equation (D6) 
can be written in the form 


_d 

dt 


=7^ < P tc,n 
c n 


’tv,n )(l + °-2M 2 ) 


-7/2 


(D7) 


Equations (D2) to (D4) and (D7) are used to represent the gas dynamics of the n th vol- 
ume. The simulation neglects the difference between static and total pressures due to 
the difficulty in computing the Mach number. Neglecting this effect results in a slight 
error in equations (D4) and (D7). Since the interstage Mach numbers within the com- 
pressor are less than 0. 5, the range of corrections for equations (D4) and (D7) are 


/ o\5/2 

1.0 < (l + 0.2M z j < 1.13 


/ 9 V 7 / 2 

0.84 < (l + 0.2M^j <1.0 

The maximum error introduced by this simplification is about 15 percent; however, this 
approximation is only reflected in the dynamic performance and does not alter the 
steady- state characteristics. This fact can be readily observed from the computation 
technique shown by the block diagram (fig. 27). The conversion from static to total val- 
ues is essentially accomplished through coefficient potentiometers 1 to 3, which are in- 
puts to the integrators. Coefficient potentiometer 6 is used to convert total to static 
density and thus merely removes the correction introduced by potentiometer 1. Since 
both integrator input potentiometers are identical, only the dynamic performance of the 
compressor will be altered. 

In the steady- state case, the expressions (D2), (D3), and (D7) become 


:/< T tc> c ,n- T tv>c,n + l> = ° 


n 


(D8) 


50 



(D9) 



( W c,n- W c,n + 1> = 0 


A n g , 

-S- (P 


- R 


tc,n tv, n 


,) = 0 


n 


(DIO) 


which implies that 


W = W i 
c,n c,n+l 


P = P 
tc,n tv,n 


T - T 
tc,n " tv, n 


(Dll) 

(D12) 

(D13) 


From equations (Dll) to (D13), it is evident that the steady- state results are not 
affected by the approximations in the volume dynamics. 

Interstage bleed can be included in the compressor simulation by modifying equa- 
tions (D2) and (D3). The bleed flow can be computed from the relation 


W b, n k b^b, n 


tv, n 


/V 


(D14) 


Including the effects of interstage bleed allows equations (D2) and (D3) to become 


“ ^sv, n T tv, J " y < T tc, n^c, n - T tv, n^c, n+1 “ T tv, n^b, n^ 


(D15) 


n 


_d 

dt 


^sv, ■ v 


( W c,n- W c,n + 1- W b,„) 


(D16) 


The inclusion of bleed flow is shown dashed in figure 27. The steady- state flow through 
the compressor is altered by the interstage bleed by forcing the flow of the (n+1)** 1 stage 
to a lower value. 

As noted in figure 27, the simplification introduced by neglecting the effects of 
static to total conversion only affects the dynamic performance of the compressor. 
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Figure 27. - Block diagram for n* h stage computation. 
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APPENDIX E 


BIVARIANT FUNCTION GENERATION 

The need arises in the compressor and turbine simulations to generate functions of 
two variables. In general, bivariant function generation is difficult with pure analog 
devices. A number of schemes varying primarily in complexity are available for analog 
use. 

The functions here required are nonlinear empirical functions of two variables, 
which are not separable into functions of single variables. We, therefore, required a 
three-dimensional surface. If the functions were separable, they could be separated 
into two single- variable functions and combined in some algebraic manner. 

Among the techniques available are correlation, curve fitting, and multicurve inter- 
polation. 

Correlation . - There is a possibility that a single- parameter correlation between 
the variables does exist. An analysis based on this premise is developed and applied to 
the turbine data in appendix F. 

Curve fitting . - Another analytical technique frequently used is curve fitting. The 
family of curves of one variable are individually fitted, and the coefficients plotted on 
diode function generators. These are then used in a circuit that generates the type of 
curve required; for example, power, exponential, and so forth. Accuracy and fre- 
quency response are very good if all electronic components can be used and interpolation 
between curves is minimized. However, a large complement of equipment is required. 

A typical coefficient plot and block diagram for a third- order power curve fit of the 
turbine enthalpy function is shown in figures 28 and 29, respectively. The curves are 



Figure 28. - Typical coefficient variation for third-order bivariant function generator. 


53 
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y = fix, z) 


Figure 29. - Block diagram of bivariant function generation using typical third-order curve fit. 


reasonably well behaved and lend themselves to easy simulation. Errors can occur 
through assuming linear interpolation between points. This could be improved by taking 
nonuniform input segments. 

Multicurve interpolation . - Diode function generators and tapped potentiometer 
servo groups can also be used to produce functions of two variables. Various interpola- 
tion schemes have been developed to increase response and reduce the amount of equip- 
ment. This method is particularly amenable to the turbine simulation because of the 
low-frequency nature of speed and the high-frequency nature of pressure ratio. 

Diode function generators are used to generate various members of the family of 
curves for constant values of the low-frequency variable. The high-frequency variable 
is their input. The simultaneous outputs of these generators are applied to the taps of 
a servo potentiometer, which is driven by the low-frequency variable. The servo mul- 
tipliers can be replaced by electronic multipliers. 

For the J85 turbine, the corrected speed is the low-frequency variable. The three 
curves for constant corrected speed were chosen to represent the enthalpy and flow per- 
formance over a range of turbine pressure ratios, the high-frequency variable. Outputs 
were applied to the end and center taps of the servo potentiometer, which was driven by 
the corrected speed parameter. 

The voltage on the taps varied inversely with the corrected speed and the servo in- 
terpolates linearly between the curves. To compensate for this, the corrected speed 
signal was modified by another function generator that alters the servo input resulting 
in a nonlinear interpolation. 

A typical block diagram for this scheme is shown in figure 30. The nonlinear cor- 
rection is given in figure 31. 
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APPENDIX F 


TURBINE PERFORMANCE CORRELATION 


A method is developed in this appendix by which the performance of a turbine can be 
analytically represented for simulation over a range of operating variables. The method 
uses the equivalence of momentum and energy as a basis. With an analytical expression 
for efficiency, the correlation is performed in terms of engineering variables and 
parameters. 

The fundamental relation of torque and rate of change of angular momentum, as 
applied to a rotating mass particle, is derived from Newton T s second law: 


L=Al Wv (Fl) 

p dtg 9 


Assuming a fixed, or average, radius of rotation and expanding the differential result in 


r / dv 0 

l p1 w 1T + v * 


dW 1 
dt } 


In difference form this becomes 


r I Av 0 

L =- W — 2 + v„ 
P g\ At 9 


AW 
At ) 


(F2) 


(F3) 


Each particle of mass passing through a point will experience this same change in 
momentum in At time. If n particles flow in time At, the resultant torque is the 
summation of the momentum change for each particle, or 


t 



nW 


Avg 

At 


+ V, 


AnW 
At t 


But, weight flow rate is defined as 


(F4) 


W = 


nW 


At 


(F5) 
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so that 


L = - (W Av^ + Vg AW) 


g 


(F6) 


For steady flow conditions 


L = — W Av fl (F7) 

g 9 

For a rotary machine, the power is given by 

& = LW (F8) 

and 

U - rW (F9) 


Combining equations (F7) to (F9) gives 

& = - W Av 


(F10) 


which is the power output in terms of angular momentum. 

Consider now the performance of a turbine from the energy viewpoint. The general 
energy equation derived from the first law of thermodynamics gives 


(P = v J — We Ti 
dt p 1 


(Fll) 


for a quasi- steady process, where the efficiency rj refers to turbine overall efficiency. 
Assuming a constant specific heat gives 


cP - 7]J 



+ C 



(F12) 


Or, as before, in difference form 


» J (k p AT t + c P T t^ 


(F13) 
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where AT; refers to the ideal total temperature drop across the turbine. This reduces 
for steady state to 


= 77<jwc at; 

p ^ 


(FI 4) 


which represents the power output in terms of energy. 

Proceeding with the isentropic relations between pressure and temperature, equa- 
tion (F14), in terms of pressure ratio, becomes 


^™ c p T t,in 


p \(r-i)/y 

J _ . t, out 

t, in ; 


(F15) 


Defining 


v s = 1/ e Jc p T t,in 

/p \(y-l)/y 

l_p,out 

r 

\ P t, in / 


as the spouting velocity, equation (F15) becomes 


W 2 


- n — v 
g 


s 


(F16) 


(F17) 


Eliminating power between the momentum and energy relations, equations (F10) and 
(F17), respectively, yields 


UAVg 



(F18) 


or, solving for efficiency 


V = 


U Ay 9 


v. 


v. 


(F19) 


Stewart, in reference 9, gives an assumed general relation of efficiency as a func- 
tion of an overall work- speed parameter: 


V = fOO 


(F20) 
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where 


u 2 _ _ /u \ 2 

gJAh 2 [vj 
s 


Hence, 



Since Av 0 /v g is not a constant, equations (F19) and (F22) imply that 



In terms of measurable variables 



(F21) 


(F22) 


(F23) 


(F24) 


(F25) 


(F26) 


This is the general functional relation for the theoretical performance of a turbine based 
on the efficiency relation of Stewart. 
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Figure 32. - Turbine torque performance correlation. 


The relation shows that the turbine performance can be expressed as a function of a 
single parameter. This makes it most amenable to simulation. 

The correlation was applied to the turbine overall performance shown in figure 11. 
The results are plotted in figure 32. The correlation is generally good. The midrange 
shows effects of a variable perhaps not included in the analysis. Some scatter can also 
be attributed to interpolation of the basic turbine data. 

Examination of the turbine flow performance shown in figure 14 suggests that the 
well known nozzle flow parameter 



(F27) 


could be used to correlate the data as a function of a single variable. 

Since a turbine consists of alternating stator blades or nozzles and rotating blades, 
expansion occurs in the stators and may also occur in the rotors. Thus, the flow 
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through a turbine stage can be represented by 



(F28) 


where P g /P^ is the stage pressure ratio. 

If the stage pressure ratio is assumed to be equal for all stages of a multistage tur- 
bine, then 



(F29) 


(F30) 


(F31) 


(F32) 


However, it has been shown experimentally that this flow correlation varies as a function 
of the speed parameter We have, finally, a correlation of the form 



(F33) 
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The results of this correlation, as applied to the J85 turbine performance data with 
n = 2, are shown in figure 33. 
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Figure 33. - Turbine flow correlation. 
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